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Abstract 

We derive a new regularization method for the calculation of the (massless) scalar self-force 
in curved spacetime. In this method, the scalar self-force is expressed in terms of the difference 
between two retarded scalar fields: the massless scalar field, and an auxiliary massive scalar field. 
This field difference combined with a certain limiting process gives the expression for the scalar 
self-force. This expression provides a new self-force calculation method. 
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I. INTRODUCTION AND SUMMARY 



The motion of an electric charged point-like object in a fixed background spacetime, is 
affected by the coupling between the object's own charge, and the field that this charge 
induces. This coupling results in a self force (also known as the "radiation reaction force") 
acting on the object. At leading order, the object's acceleration due to this self force (in 
the absence of non-gravitational external interactions) is proportional to q 2 / fi, where q 
and fj, denotes the object's charge and mass, respectively. This leading order is obtained 
by treating the particle's field as a linear perturbation over a fixed curved background 
spacetime. Analogous to the electromagnetic self-force there are other types of self forces: 
a scalar self-force is induced by a scalar charge, and a gravitational self-force is induced by 
the object's mass (in this case the object's acceleration due to the gravitational self-force is 
proportional to //). 

In recent years, there has been growing interest in the calculation of the self force in curved 
spacer t fl Form a. ^. of the 8 e!f foree have been derived: Mine, Sasa kl and 
Tanaka [1], and independently Quinn and Wald 0, have rece ntly obtained an expression for 
the gravitational self force; previously, DeWitt and Brehme jl5j obtained an expression for 
the electromagnetic self force; and Quinn recently obtained an expression for the scalar self 
force. In the case of a weak gravitational field, explicit self-force expressions were found by 
DeWitt and DeWitt Q| , and by Pfenning an Poisson 4] . Analysis of the self force in curved 
spacetime also has a practical motivation: one possible source for LISA - the planned space- 



based gravitational wave detector [17|, is a binary system with an extreme mass ratio, which 
inspirals toward coalescence. Here, the self force is required for the calculation of the accurate 
orbital evolution of such systems. These orbits are needed in order to design templates for 
the gravitational waveforms of the emitted gravitational radiation. A calculation method 
for the self force in such binary systems, was devised by Barack and Ori 0, this method 
was recently improved f], 0, and also implemented numerically in certain cases 0,0. For 
other approaches to the self force problem see [3, [ll|, [l^ . 

Previous analyses have provided several expressions for the scalar self force, which are 
equivalent to each other. In Quinn's derivation 0, the expression for the scalar self force 
is composed of two different types of terms: explicit local terms; and a non-local term, 
expressed as a certain integral over (the partial derivatives of) the retarded scalar Green's 
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function. We shall refer to these terms as the local terms, and the non-local term, respec- 
tively. Very recently, Detweiler and Whiting [14|, developed a different method to express 
the self force in curved spacetime. In the scalar-field variant of their analysis, they showed 
that the scalar self force can be obtained from ip R , which is a certain non-retarded solution 
of the homogeneous scalar-field equation. 

In calculating the self force acting on a point like object, one typically encounters a 
divergent expression, from which the finite (and correct) expression for the self force should 
be obtained using a certain regularization method. In this manuscript we present a new 
self-force regularization method. We consider a point-like scalar charge in curved spacetime, 
which induces a scalar field 0. We show (under certain assumptions) that the scalar self 
force acting on this scalar charge can be expressed in terms of two retarded scalar fields, 
which satisfy prescribed differential equations. These fields are: <fi - which satisfies the 
inhomogenous massless scalar field equation, with a charge density p; and m - which satisfies 
the inhomogenous massive scalar field equation, with exactly the same charge density p (here 
m denotes the mass of the massive scalar field). More specifically, we show that the scalar 
self force is given by 



Here, zq is the self-force evaluation point on the object's world line, and x is a point near 
the word line, defined as follows. At z we construct a unit spatial vector n M , which is 
perpendicular to the object's world line but is otherwise arbitrary (i.e. at z Q we have n^n^ = 
1 , n^u^ = 0). In the direction of this vector we construct a geodesic, which extends out 
an invariant length 5 to the point x(zq, n^, 5); throughout this manuscript and denote 
the object's four-velocity and four-acceleration, at Zo, respectively. 

We show below that the value of f^ el ^ given by Eq. {1} is identical to the scalar self force 
obtained by Quinn 3|. Therefore, the problem of calculating the scalar self force in curved 
spacetime is equivalent to the problem of solving the two scalar partial differential equations 
for the fields and (f) m , and then carrying out the prescribed set of operations summarized 
by Eq. (0). Note that for the massive field <f> m , only the asymptotic behavior as m — > oo is 
required. 




(1) 



where, 




(2) 
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We mention here a similar method developed by Coleman [18[ for the derivation of the 



□ 



electromagnetic self force in flat spacetime (the Abraham-Lorentz-Dirac term |l9l l20l. l2l|). 
In this method, the self-force regularization is obtained by replacing the Green's function 
of the electromagnetic four-potential with a different (more regular) function, which de- 
pends on a certain parameter. This new function behaves like the original electromagnetic 
Green's function at the limit where the parameter approaches infinity. Coleman obtained 
this new function from the Green's function of a fictitious massive electromagnetic field. 
This regularization method is similar to Pauli-Villars regularization in quantum field theory 

Q 

The derivation of Eq. (JIJ is based on several properties of the fields and m , which we 
now briefly summarize. Previous methods 0,0,0,13 used the following decomposition of 



the scalar field near the object's world line: 0(x) = <fidir(x) + ^>taii{x)- In this decomposition 
the direct field 4>dir(%) is a field that propagates along the future light cone emanating from 
the object. Near the object's world line this direct field can be determined using a local 
analysis jl^ . and it diverges as 5 — ► 0. The second field in this decomposition, the tail field 
4>taii{x), is of nonlocal nature, i.e., its value generically depends on the entire past history of 
the object's world line. Moreover, the tail field is finite at the limit 5 — > 0, and it vanishes in 
flat spacetime (in 1+3 dimensions). This tail field was found to be useful in the calculation 
of the self force in curved spacetime p - the above mentioned nonlocal term of the self 
force, can be directly obtained from (fitaii,^ near the object's world line. 

In a similar manner, near the object's world line the massive field <f> m (x) can be decom- 
posed according to <f> m (x) = 4> m (dir)(x) + (prn(taii)(x). Here again, the massive direct field 
4>m(dir) is a field that is propagated along the future light cone emanating from the object, 
and it diverges as 5 — > 0. The massive tail field <fi m (taii), is a field of a non-local nature, and 
is non-singular on the world line (for a finite value of m). However, unlike the massless tail 
field, this massive tail field does not vanish in flat spacetime. 

The formal expressions for the two direct fields <fi m (dir) and <pdir (see below) show that they 
are identical. This property is an important element in our construction, since it implies 
that by subtracting <p rn from <p [see Eq. (j2J)] we obtain 

A(f)(x) = (j) ta il{x) - <f>m(tail){x) , (3) 

which is finite as 5 — > 0, while the individual fields <p and (f) m diverge at this limit. Similarly, 
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the partial derivatives A0 iAt are also finite at this limit. However, their values at this limit 
depend on the direction along which the limit 5 — > is taken (i.e., they depend on n^). 
This directional dependency is completely annihilated by the term (qm 2 /2)n fl in Eq. (Q). 
Therefore, the expression obtained by taking the limit 5 — > in this equation does not 
depend on n^. To illustrate some of these properties consider the following simplest example. 
Consider a static unit point charge situated at Zo, in a flat spacetime. In this case, the fields 
4> and 4> m , which satisfy equations (jU and (JTj), respectively, are given by 

= 4>dir = ~ , 4>m = ~ + -(-1 + e~ mr ) . 

Here r = \x — Zq\. We therefore find that A<f> = —4> m (taii) — r(l — e~ mr ), and the gradient of 

2 

this expression in the vicinity of the particle is V(A0) = —^r + 0(r), where f is a radial 
unit vector. Hence, both A0 and its gradient remain finite at the limit 5 — > 0. Moreover, in 
this case the self force obviously vanishes - a result that is easily derived from Eq. (|T]). 

The calculation of the scalar self force, naturally involves derivatives of the field - A0 
in our construction. By virtue of (J3j), these derivatives will include both 4>tau,fi an d (pmitaii),^- 
As we mentioned above, the non-local part of the self force can be obtained from (/)t a u,n- 
Therefore in our method 4> m (t a ii),^, should give rise to the (above mentioned) local terms of 
the self force 1 . This might look strange at first sight, because the derivatives 4> m (taU),fi depend 
on the entire past history of the object world line. Note, however that in our method the 
self force is obtained by taking the limit m —>■ oo in Eq. ((TJ. At this limit the contribution 
to Eq. (|TJ) from the term 4> m {taii),^ is a local one (see calculations below). This behavior may 
be traced to the properties of the massive scalar Green's function at the limit m —>■ oo. 

In Sec. iniwe present the massive-field approach, and show that f^ el ^ which isgiven by 
Eq. (JTJ) is equivalent to the standard expression for the scalar self force, given in [3[. 

II. THE MASSIVE FIELD APPROACH 

We consider a scalar field <p(x), which satisfies the scalar field equation 

□0 = -4tt p . (4) 

1 We comment that there is also a certain local contribution to the self-force coming from 4>taii,^', this local 
contribution is taken into account in the calculations below. 
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Here D(p = and p(x) is the charge density of the scalar object. We use the signature 
(-+++), and natural units where c = 1 throughout. Our purpose is to derive an expression 
for the self force on a point-like object. We therefore consider a charge density of a point 
particle. The particle's world line is denoted by z(r), where r is the particle's proper time 
(we allow an arbitrary acceleration - presumably caused by an arbitrary external force acting 
on the particle), the particle's scalar charge is denoted by q and the particle's charge density 
is therefore 

f°° 1 

p(x) = q / 5 4 [x — z(r)]dT . (5) 

J-oc v-g 

Here g(x) is the determinant of the background metric. Next, we define the field of the 
scalar force to be 2 

F^x) = q<p^ . (6) 

This field diverges as x approaches the particle's world line. In order to extract the self 
force from this singular field, we need to regularaize this expression. For this purpose, we 
introduce an auxiliary field <f) m (x) which satisfies the massive field (Klein-Gordon) equation 3 

(□ - m 2 )0 m = -Aup . (7) 

Here the charge density p is the same charge density introduced above [i.e., it is given by 
Eq. @]. First we consider the field's mass m to be some fixed (large) quantity. Similar to 
Eq. © we define the field of the massive scalar force to be 

F(m)p(x) = # m>/r (8) 

Like the massless field the massive field F^^ diverg approaches the particle's 

world line. However, as we show below, their difference remains finite at this limit. We 
therefore construct the difference field A<j) given by Eq. Q, and introduce the difference 
between the massless and the massive scalar forces 

AF^(x) = F M (x) - F (m> (x) = qA^ . (9) 

This field is an essential element in our construction, since in our method the scalar self 
force is calculated from AF^ - see Eq. (0). To analyze the properties of this field, we derive 



2 This definition conforms with definitions given in a different definition is also possible ^3 . 

3 Note that here the field's mass m has the dimensions of (length)^ 1 . 
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a formal expression for its value near the particle's world line. This expression will later be 
used to obtain Eq. (JTJ). 

First, we derive a formal expression for the field A(p using the corresponding retarded 
Green functions of the massive, and massless fields. The retarded solutions of equations 
(J4l7j) for the charge density given by Eq. (j5J), are formally given by 



cj)(x) = q G(x\z(r))dr, (10) 

J —oo 

POO 

4>m{x) = q / G m (x\z(r))dr . (11) 



Here, G(x\x') and G m (x\x') are the corresponding retarded Green's functions which satisfy 

4-TT 

UG{x\x') = ~-^=5\x-x'), 
Ait 

(□ - m 2 )G m {x\x') = --^=5\x-x'). 

Both Green's functions vanish for all points x which are outside the future light cone of x'. 
We assume that the spacetime is globally hyperbolic, and therefore G and G m are uniquely 
determined by the above requirements. We further assume that the integrals in equations 
(|1UI11|) converge in some local neighborhood of the particle's world line, but not on the 
world line itself. Expressing A<p with these Green functions gives 

/CO 
AG{x\z{r))dT , (12) 
-co 

where AG(x\x') = G(x\x') — G m (x\x'). Consider the above expression for A<ft(x) in some 
local neighborhood of the particle's world line. In the discussion below we will use two 
separate expressions for AG(x\z) in Eq. ()12j) : a local expression - for the case where z is in 
a local neighborhood (defined below) of x, and a non-local expression for the case where z 
is outside this local neighborhood. In the next subsection we derive the local expression for 
AG. 

A. Local expression for AG 

In the vicinity of any point x' there is a local neighborhood, in which any two points can 
be connected by a unique geodesic within this neighborhood, see theorem 1.2.2 in (this 
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neighborhood is sometimes called a geodesically convex domain) . In this neighborhood, the 



massless and massive retarded Green's functions are given by 



24] 



G(x\x) = 0(E(x),x')[U(x\x')5((j) - V(x\x')e(-a)\, (13) 
G m (x\x') = 0(E(x),x')[U(x\x')5((j) - V m (x\x')G(-a)\. (14) 

Here, a = a(x\x') is half the square of the invariant distance measured along a geodesic con- 
necting x and x'] a is negative for a timelike geodesic, positive for a spacelike geodesic, and 



vanishes for a null geodesic; U 
definitions and properties, see 



x\x'), V(x\x'), and V m (x\x') are certain bi-scalars (for their 
23]); S(x) is an arbitrary space-like hypersurface containing 
x; and O (£(#), x') equals unity if x' is in the past of and vanishes otherwise. 
From equations (|13ll4j) we find that locally (i.e., in a geodesically convex domain) 



AG(x\x') = -AV(x\x')Q(J:(x),x')e(-a) , (15) 

where, 

AV(x\x) = V(x\x') - V m (x\x) . (16) 

Note that the same "direct term" U(x\x')S(a) appears in both Green's functions in equations 
(|T3*|) and ([14)1. and therefore it cancels upon their subtraction. 

Each of the direct fields 4>dir and 4> m (dir), which were mentioned in the previous section, 
is obtained by integrating over the corresponding direct term. These fields are therefore the 
same, and are given by 



4>dir{x) = 4> m (dir){x) = q 0{Y,{x),z)U(x\z)5{a)dT . (17) 

J T~—E 

Here r~ denotes the retarded proper time [i.e., the proper time at the point of intersection 
of the past null cone of the field evaluation point x, with the particle's world line z(r)], and 
e is an arbitrary small time interval. Since these direct fields are equal, they cancel each 
other in Eq. (J2J), which gives Eq. (J5J). 

It will be useful later to have the explicit dependence of V m on m; for this purpose we 
use the Hadamard expansion. It was shown by Hadamard [2^ that the bi-scalars V and V m 



S 



can be expanded as follows 4 : 



V(x\x') 



—rv n (x\x') , 



(18) 



n=0 



V m (x\x') 



y2—rV„(x\x') , 

' n! 

n=0 



(19) 



where the coefficients f n and satisfy certain differential equations. In the case of an 
analytic metric, for every point x' there is a local neighborhood in which the Hadamard 
expansion is uniformly convergent [2£| (This is also true for a class of non-analytic metrics 
|26j.). The local neighborhood of x' in which both equations (|18ll9j) converge uniformly 
is denoted by D(x'). Here and below we shall consider the Hadamard expansion inside 
D(x f ). Expressing the coefficients v n in terms of the coefficients v n , and the bi-scalar U 
gives [The derivation of this expression is given in Appendix A, for a different derivation, 
see Friedlander [23] equation (6.4.19).] 



Here J n denotes the Bessel function, and we introduced s = \J —2o . Note that the bi-scalar 
V m is an even function of s. We point out that the local expression for AG inside D(x') is 
obtained by substituting equations 1)18120)) in equations ()15I16|) . 

Before continuing the detailed calculation, we briefly discuss some of the properties of V m , 
and their relations to the self force in our method. In the prescription given by Eq. (Q) the 
asymptotic expression of lim^o A0 jAt , as m — * oo is required. Now, Eq. ()12J) implies that 
A0 can be determined by an integral over AG, and by equations ()15ll6j) the local expression 
of AG depends on V and V m . Therefore, to obtain the asymptotic form (as m — > oo) in 
Eq. we have to study the asymptotic form of V m (as well as the asymptotic form of the 
nonlocal expression of AG). From Eq. (|20j) we find that as m —>■ oo the Bessel functions 
oscillate rapidly (for s ^ 0); these rapid oscillations have a canceling effect upon integration 
with respect to s. In fact, the calculations below show that at the limit m^ooan integral 
(with respect to s) over the first term in Eq. (j2"Uj) behaves much like an integral over the 
direct term U5(cr). This is an important property in our construction 5 , because it explains 

These two expansions are defined slightly different from the various definitions in 0, 0, and 0- The 
coefficients of these various definitions are related to the coefficients in equations (|f 811 9fl by multiplying 
by certain numbers. 

5 A similar property was used by Coleman in the calculation in flat spacetime |l8| . 




(20) 
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the source of the local terms of the self force in our method. In other methods, for example 
in , the local terms of the self- force are essentially obtained from the value of (pdir ^ near 
the world line, by some regularization method. Here however, in the expression for A0 the 
direct fields (pair and (f> m (dir) canceled upon subtraction. Therefore it might seem strange 
that we still obtain the local terms of the self-force. This is possible because at the limit 
m — > oo we find that the first term in Eq. (J2*U|) behaves similarly to the above mentioned 
direct term. 



B. Calculation of AF^ near the particle's world line 

We now use the above local expression for AG, and derive an expression for the field 
AF tl (s) near the world line z(t). Before proceeding, we define some notation. The particle's 
proper time at the self force evaluation point z is denoted by r ; The past intersection of 
the boundary of D(zq) with the world line is denoted by z{t\). By our construction the 
points z and any point z(f), such that ro > f > Ti, can be connected by a unique geodesic 
within D(zo). It may be that all these geodesies are timelike geodesies . If this requirement 
is not initially satisfied, we then increase the value of T\ such that all these geodesies with 
tq > f > T\ will become timelike geodesies. The field evaluation point x is taken to be 
sufficiently close to the particle's world line (i.e. 5 is sufficiently small), such that both 
points x and z(t~) are within D(zo) and r~ > n. 

Using equations (|12|15|h we obtain 

A(j)(x) = -q AV(x\z)dr + q AG(x\z)dr. (21) 

J T\ J — OO 

This field has a finite value as x approaches the world line (see Appendix[nj). From equations 
flUEO we find that 

AF^x) = -g 2 (r-), M [AVV " <f f ^Vdr + q 2 T AG /r . (22) 

Throughout this manuscript the indices fi, v refer to the field evaluation point (here this 
point is denoted by x), and the subscript r~ indicates that the variable z(t) inside the 
brackets is evaluated at the retardation point z(r~~). Using equations (jl6pi8|2()}) we obtain 

[AV(x\z)] T - = -\^ 2 [U(x\z)] T - . 
10 



This relation together with Eq. (|22J1 give the following expression for AF^ 

AF^(x) = F taU ^(x) + 1 -q 2 m 2 [U] T -(r-) >fi + q 2 f V m ^dr - q 2 P G m Jr. (23) 

* Jti J —oo 

Here we introduced, 

F tail ^x) = q 2 ! 1 G^dr -q 2 I V^dr. (24) 

J —oo Jt\ 

It is useful for the calculations below, to change the splitting of the integration interval 
in Eq. into a "smooth splitting" defined as follows. We introduce auxiliary weight 

functions h\{r) and h2{r) (sufficiently smooth, see below), and an arbitrarily small time 
interval e, such that h\{r) = 1 for r > T\ + e, hi(r) = for r < T%, and hi(r) varies 
smoothly in the interval [ti, ri + e]. Defining /^(t) = 1 — hi(r), we find from Eq. (}2*3*j) that 

AF^(x) = F ta \(x) + \q 2 m 2 [U] T -(r~)^ + q 2 [ h x V m ^dr - q 2 [ ' + h 2 G m ^dr . (25) 

^ Jri J -oo 

Eq. (J25)) provides a general formal expression for the field AF ^ near the object's world line. 
Next, we shall use Eq. (J25J) to show that Eq. is identical to the scalar self force. 

C. Derivation of the self force expression 

We now follow the prescription given by Eq. (P), and perform the following successive 
operations on the field AF^: 

i Calculating the limit 5 — > of AF M . 

ii Calculating the asymptotic form of lim^o AF ^ as m approaches infinity. 

These mathematical operations will be performed separately on each term in equation (|25|). 



1. The first term 

Consider first the limit 5 —>■ of the first term F tad M in Eq. (|23|) . Since V is smooth (see 



theorem 4.5.1 in 



23j), and r (x) — > r as 5 — > 0, we find from Eq. (|2*4^) that 
n /-to 



/Tl /-To 
GJz \z(r))dr-q 2 j V^z Q \z{r))dr . (26) 

This expression is equivalent to the non-local term of the scalar self-force which was found 
in Ref. jjj . Since, this term is independent of m, it is not affected by limit m — > oo. 
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2. The second term 

Consider next the second term in Eq. (|25|). The second term can be expressed as 6 

\q 2 m 2 {U) T - (r-)^ = -\<?m\U) T - . (27) 

Here the index a refers to the point z a (r) on the particle's world line, and u a = The 
calculation of the limit 5 —>■ of this expression follows from a local expansion of the various 
terms in Eq. (j2Zj). The detailed calculation is given in Appendix B, there we show that 

\im^q 2 m 2 (U) T -(T-) >fi = ~q 2 m 2 (u M + n M ) . (28) 



3. The third term 



Consider next the third term in Eq. (|25)l . First, we take the limit 5 — > of this term. 
Similar to Eq. (pjj) we find that 



TO 



lim / hiV^ „(x|z(r))dr = / h x V m Jz \z(r))dr . 

We now split the integration interval into two intervals: one interval is \t\ + e, r ] in which 
/?,! = 1, and the the other interval is [n, n + e] in which varies smoothly, this gives 

TO 



hV m ^(z \z(r))dr = V m ^(z \z(r))dr +/ /i i V r m ^(2;o|z(r))rfr . (29) 

First we focus on the first integral on the right hand side of Eq. (|2*H|) . We change the 
integration variable to s (with z fixed). Recall that for any point z(f) on the world line 
such that n < t < To, we have a unique timelike geodesic within D(zq), which connects zq 
and z(f). For each 2;(f) in this range we have 

^(z \z(f)) = s, a ^<0. 

This expression is negative because s jCf is a timelike (future directed) vector at z(f) which is 
tangent to the geodesic that connects z(f) and z ; and ^ is the four velocity vector (also 
future directed) at z{f). 



6 This expression follows from the relation between r and x on the lightcone, which is given by the equation 
ct(x\z(t~)) = 0; see for example 
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Substituting Eq. (|2U|l into the first integral in Eq. 
integration and summation gives 

Vm,Jzo\z(r))dT = 



and interchanging the order of 



ri+e 
oo 



£ 

n=-l 



— m) 



v n ^J n {ms)s n + [J„(ms)s n ] 



dr 
ds 



-ds 



Here we introduced s\ = s{zq\z{t\ + e)), u_i = [/, and the various bi-quantities are evaluated 
at the points z and z(r) [e.g., s = s(2b|z(T)), v n = v n (z Q \z(r)) etc ...], this notation is used 
in the rest of this subsection. Using integration by parts we find that 



V m ^{zo\z{r))dT 



(30) 



ri+e 
oo 



^ / Jnijns) 



S 



m' 



dr 
ds 



V n ... - S, 



* ds 



d (dr 



ds + 



Here, in the differentiation with respect to s, Zq is fixed; whereas in the differentiation with 
respect to x M , z a (r) is fixed. Local analysis shows that the first boundary term in Eq. (|30|) 
is given by (see Appendix B) 

' m 2 R 
~2~ 



Vm Ts s ^ 



s=0 



. u,, . 
12 ' 



(31) 



Here R denotes the Ricci scalar 7 . We do not need to calculate the second boundary term in 
Eq. (J3U|) since it cancels with another term (see below). 

Next, we calculate the asymptotic form of the integral in Eq. ()30|) as m — > oo. Consider 
the integral over the n'th term in this equation. To calculate the asymptotic form of this 
integral, we expand the term inside the square brackets in powers of s. Integration over the 
terms in this expansion requires one to calculate integrals of the following form: 



o 



JJms) 



-ds 



77V' 



2n+Q+l 



J n (y)y n+a dy. 



(32) 



msi 



Note that the coefficients v n in Eq. (|3~U|l are all smooth functions (see 23], sec. 4.3), and 
the expansion of the expressions that do not depend on v n , contain only positive powers of 
s (see Appendix B). Therefore, in Eq. ()32|) we need to consider only a > 0. As m — ► oo 



7 Here the Riemann tensor is defined with the opposite sign with respect to the definition in |lo| . 
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the integral in this equation vanishes for n > 0. Therefore, only the terms with n = — 1 
contribute to Eq. (j3*Uj) . In this case (n = —1) the expansion of the term in the square 
brackets in Eq. (J3U)) includes only integer powers of s (see below). From these terms only 
the terms with a = 0, 1 do not vanish as m — > oo. Therefore, we need to expand the n = —1 
terms only up to the first order in s. These expansions give (see Appendix B) 

TT d ( dr\ 1 1,. 2 \ ^/ 2\ /oo\ 
ds V ' M ds J = ~2 apL + 3^~ a u v) s + °( s )> ( 33 ) 

^ - s,J^ = -~ {R/u v + itT W ,) s + 0(s 2 ) . (34) 

Here d u denotes the covariant derivative of a u with respect to r. Here and below (unless ex- 
plicitly indicated otherwise) the coefficients u^, a u , d^, W v ', R are evaluated at z . Equations 
fl3()l3H32l33l34l) give 

V m Jz \z(r))dr^ (35) 



ri+e 



^m 2 u u - ^ma u + -(d M - a 2 u^) + (-R^ u u v + -RP u u v u v u u - 



V -s 

^ / 

Here, and throughout this manuscript, the symbol = represents equality up to terms that 
vanish as m — > oo. 

We now focus on the second integral on the right hand side of Eq. f!29j) . Using a 
calculation similar to the one that was performed on the first integral in Eq. ()29j) . we find 
that 

hiV m ^(z \z{T))dT = (36) 



n=-l" /Sl L v 



dr, / dv n \ d ( dr 

h\ v n — s — — — — hi—s 

\ 'ha ,a Q s ds\ d.S 



ds + 



hlVm Ts'^)^ -{ Vmhl Ts s >» 



=si 



Here we introduced si = s(zo\z(tx)). Recall that hi(s\) = 0, and hi(§i) = 1. Therefore, 
the second boundary term in Eq. (JBljj) vanishes, and the first boundary term is minus the 
boundary term in Eq. (J33j) . Noting that asymptotically J n {x) = J ^cos[x — (n + |)|] + 



0(x~ 3 / 2 ) |23, we find that as m — > oo the integral in Eq. (|36|) vanishes for n > 0. For 
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n = — 1 we substitute J_i(ms) = —4-Jo(ms) in this integral, and integrate by parts once 



more. Taking the limit m->oowe find that s 



hiV m ^(z \z(r))di 



V -s 

m ds 



(37) 



From equations 



we obtain 



hiV m Jzo\z(r))di 



(38) 



+ {\R^u u + l -R^u v 





12 



Koto that the local t(Tnis of the scalar self force found in Ref. j^j arc identical to tlic terms 
that are independent of m in Eq. 



4- The fourth term 

Consider next the fourth term on the right hand side in Eq. (|23|) . Taking the limit 5^0 
of this term gives 

"Ti+e 

h 2 G m ^{z \z{T))dr . (39) 

We now discuss the limit m — > oo of this expression. Notice that in the above calculation 
of the third term, we found that the value of the third term at the limit m — > oo does not 
depend on T\ [see Eq. f!38|) ]. This suggests that the entire fourth term should vanish at the 
limit of interest. To further support this statement we introduce the following physically 
motivated assumption. It is well known that in flat spacetime as m — > oo the range of the 
massive field interaction approaches zero - here we shall assume that this statement is also 
valid in curved spacetime. Therefore, the massive scalar force F{ m )n f° r an y point that is 
not on the world line vanishes as m — > oo. We now consider a different world line z(f), 
which coincides with z(t) for r < t% + e, and is slightly displaced with respect to z{r) for 
r > T\ + e, such that z is not on z{j). The massive force field F^^zq) for the world line 
z is given by 

F( m )n(zo) = q 2 d IJ , / G m (z \z(f))df 

J — oo 



Here we require that hi(s) is a C 2 function. 
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where r is the retarded proper time on z. Using the same smooth splitting as in Eq. (|25jh 
we find that 

h 2 G milx dr + q 2 d fl / h x G m df . (40) 

-OO J T\ 

Notice that the first term in this equation is identical to Eq. (J39j) . By calculating the limit 
ttl — ► oo of the second integral in Eq. (J4*U|) one finds that this term vanishes at this limit 9 . 
Therefore, Eq. ()40|) together with our above assumption imply that the entire forth term 
vanishes at the limit of interest. We verify this conclusion for the special case of a de Sitter 
background spacetime (see Appendix D). 



5. Result 



By substituting the above four terms in equation (|23|) we obtain 



hmAF II = F ta \(z )+ (41) 
q 2 { -\m 2 n^ - \ma^ + ^(d^ - a 2 u^) 



The term in the curly brackets in Eq. (|4*Tj) has a directional dependence on and it 
diverges as m — > oo. To remove this directional dependence, and the divergence behavior as 
m — > oo; we subtract the term q 2 (—^m 2 n^ — -jma^) from both sides of this equation. We 
can now safely take the limit m — > oo, from which we obtain Eq. 

/f ^ Oo) = q lim { lim A<^ + \qm 2 n^ + \qma^ \ = (42) 

^ m— >oo I <5— >0 2 '2 J 



-(a^ - a 2 u^) + y^RpUv + ^R vu u v u u u^ - j^Ru^ 



This last expression is identical to the scalar self force expression found in 



13] 
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APPENDIX A: COEFFICIENTS RELATIONS 



Here we express the coefficients v n (x\x') in terms of the coefficients v n (x\x') and the 
bi-scalar U(x\x r ). From these relations we then derive Eq. (J2(Jj) . 

We consider x to be within D(x'). The coefficients v n (x\x') and v n (x\x') satisfy the 
bllowing recurrence differential equations along the geodesies which connects x and x' (see 
2J, Sec. 4.3): 



^-^) <y? = -\(n-m 2 )u. 



1 



v n (n + 1) + v n:ll - —U-^ o.v = --(□- m )v n -\ . 



(Al) 
(A2) 



U~ ift J' ' 2 

Here the index [i refer to derivatives with respect to x M . The differential equations for the 
coefficients v n (x\x') are obtained by substituting m = in equations (|Al|A2j) . We introduce 
an affine parameter r along the geodesic x(r) extending from x' to x, where x(0) = x' and 



x(l) = x. By substituting o.j 1 
integrated j23j], which gives 

v (x(r)\x') = 
v n (x(r)\x') = 



dr ' 



into equations ()A1IA2|) . these equations can be 



U 
2r 



{U-m 2 )U 



U 



U 



2r n+1 



dr' 

x(r') 

□ - m 2 )v n „i 
U 



dr'. 



(A3) 
(A4) 



x(r') 



Similar expressions for v n (x(r)\x r ) are obtained by substituting m = in these equations 
which gives 



v (x(r)\x') 



v n [x\r x 



U 
2r 



nu 
IT 



u 



2r n+l 



in 



dr' 

x(r') 
OVn-1 

u 



dr'. 



x(r') 



First, we consider vq. From equations (jA3IA5|) we find that 

v = v + ^Um 2 . 

Next we consider v\. From equation (|A5|) we find that 

1/2UU-VQ d (v r 2 



(A5) 
(A6) 

(A7) 
(A8) 
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By substituting Eq. (|A7|) into Eq. (jA4|) with n = 1 and using Eq. (|A8|) . we obtain 

v o 2 U A 

V\ = v\ + —m + — m . (A9) 
2 a 



By repeating this process of substitution of the coefficients Vi in Eq. (jA4j) for Vi+i, we find 
that v n is given by 



n / o \ n—k 



■=E(t (^f. (A«» 



fe=-l 

where w_i = [/. The Hadamard expansion for in Eq. ()19|) now takes the form 

oo „ n / 9 \ n— k 

n! 

n=0 fc=-l 



We introduce p = n — k, and rearrange the order of the terms in Eq. (jAll|) 10 which gives 
= f; A»f] ( ' 1 , ) ; ( '"ff - = mU J -i^l + ± vMms) (=£)". (A12) 



Here s = \/—2a. 



APPENDIX B: LOCAL ANALYSIS 

Here we derive equations ()33l34l31l28j) . First, consider the derivation of equation (J33j) . 
This equation requires a local expansion of the expression U 4-{s^^-) in powers of s. We 
start by expanding <j-^{zq\z{t)) for a point z(t) which is in the vicinity of point Zq, this gives 

a,,(z \z(r - t)) = [a J - t[d T a.J + £[%*.J - j[fy.J + 0(t% (Bl) 

Here we introduced the proper time difference t = r — r, and the following notation for 
the coincidence limit [er ;/i ] = lim^o c^- Here and below the indices /x, i/, 77 refer to the point 



= Zq, an d t 



he index a refers to the point z a (ro — t). Employing the coincidence limits 
given in 28 1 we obtain 

t 2 t 3 

a.^(z \z(r)) = t Ufl - -ap + -d„ + 0(t 4 ) . (B2) 



10 We do not discuss the question whether it is permissible to rearrange the order of the terms in this sum. 
We comment however that our resultant expression is equivalent to an expression for V m obtained by 



Friedlander in a completely different manner, see [2^ equation (6.4.19). 
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Here and below, unless indicated otherwise the coefficients u^, a M , a M are evaluated at zq. By 
differentiating the normalization u At (r)u At (r) = — f twice, we find that a 2 = —a^u^. From 
this relation together with s 2 = —a.^a.^ 11 and Eq. (jB2|) we obtain 

s(z \z(r))=t + ^a 2 t 3 + O(t'). (B3) 

From this equation we find that 

t = S -laV + 0( S 4 ), (B4) 

and therefore 

! = - 1 + I V + (A (B5) 
From equations (jB2|B3|B4j) we obtain 

s tfM (z \z(r)) = s-% = -u, + S -a, -s 2 (^-^+ 0(s 3 ) . (B6) 

A local expansion of the bi-scalar U gives 

U = 1 + ^R^(z )a.^ + 0(s 3 ) . (B7) 

From equations ()B5IB6IB7|) we obtain Eq. (jBlIj) 

d dT \ s 
U d's i ' S ^ = ~2 a>M + 3 ( ^ ~ ^ + ° (s2) • (B8) 



Next, consider the derivation of equation (j34j) . Here our purpose is to expand the expres- 
sion ^ [y,n ~ s ,[il^j ^ n P owers °f s - Recall that in the differentiation with respect to s, 
is fixed; whereas in the differentiation with respect to x M , z a (r) is fixed. Using equations 
()B5IB6IB7|) . and the above mentioned coincidence limits we obtain Eq. 
dr / dU\ 1 



ds 



- = ~l ( + R^u^) s + 0(s 2 ) 



where R vu is evaluated at Zq. Next, consider the derivation of Eq. (|3ip. From Eq. (|20|) we 
find that 



9 

m 



(V m ) s=0 = —(U) s=Q + (v ) s=Q . 



Employing Eq. (|B7|) . and noting that (t> ) s =o = —R/12 151 ]. we obtain 

m 2 i? 



(Kn)s=0 



12 



11 For this relation see Note, that here we defined s = \/—2a. 
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where R is evaluated at zq . From this equation together with equations (jB5IB6j) we obtain 
Eq. (jSU) 



V -s 
m ds 



s=0 



m 

T 



R 



Next, consider the derivation of equation This equation requires the calculation 

of the limit 5 — > of the expression (t~) = — [a- !fJl (a- a u a )~ ] t _. We therefore expand the 
expression inside the brackets in powers of 5, keeping only the first leading order in this 
expansion. Consider first the term a ;IM (x\z(r)). We choose the points x(t , n M , 5) and z(t) 
to be in the vicinity of the point z(tq), such that S is of the same order as t. Expanding a ;tl 
we find that 



a^(x(T ,n^,5)\z(T -t)) 



dz a 



+ 5 



dx v 

~dJ a ^ v 



0{8 2 ) 



(B9) 



where the square brackets denote the coincidence limit 



[a.^] = lim lim o.^ . 



Taking r to be the retarded proper time r = r gives (see |l 



t = 5 + 0{5 2 ) , {a, a u c 



5 + 0(S 2 ). 



(BIO) 



From equations (j27IB7IB9IB10|) together with the above mentioned coincidence limits we 
obtain Eq. 

I i 

lim-g 2 m 2 [f/] r -(r") iAt = --q 2 m 2 ( Ufl + n M ) . 



APPENDIX C: REGULARITY OF A<j>{x) AS 5 -» 



Here we show (under certain assumptions; see below) that A0(x) remains finite as 5 — > 0. 
We consider a finite value of to; hence the bi-scalar AV = I 7 — V^, is a smooth function (see 
theorem 4.5.1 in 2^|). Therefore, the first integral in Eq. (|21|) is a continuous function of S 
at z . 

Consider next, the limit 5 — > of the second integral in this equation. The non singularity 
of this limit follows from the following argument. First, examine lim^o A0(x) for a particle 
with a different world line denoted by z'(t'), and defined such that it coincides with the 
original world line z{r) for r < ti, and is slightly displaced with respect to z{r) for r > T\, 
such that is n °t on ^'^ As before, we assume that for the world line z'(t') the integrals in 
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Eq. (jlOlllj) converge in some local neighborhood of the z'(t'), but not on z'(t') itself. Since 
Zo is not on z'(r'), we find that for the world line z'(r'), A<P(zq) is finite. This implies that 
by expressing A<f)(z ) [or lirn<5_>o A<f>(x)] by Eq. (|21j) ; the second integral in this equation 
is finite for the world line z'(t'), and therefore it is also finite for the original world line 
z(t). We therefore conclude that A(p(x) is finite at the limit 5 — > 0. Similarly, one can show 
that AF^x) is finite at the limit 5 — > 0; this limit, however, has directional irregularity (it 
depends on This property is discussed in section ITT1 



APPENDIX D: CALCULATION OF THE FOURTH TERM IN DE SITTER 
SPACETIME 



Here, we provide an example of the calculation of the fourth term. We calculate the limit 
m —>■ oo of Eq. (J39)) in de Sitter spacetime. In this case the retarded Green's function of a 
massive scalar field is given by |29| 



G m (x\x) = Q(E(x),x') 



As 



sinh(As) 



3/2 

5(a) + A 2 P^[cosh(As)]0(-a) \ . (Dl) 



Here a = a(x'\x), the Riemann tensor is given by R^ UV k = A 2 (g w (/v K — g^Quri), PL is the 
derivative of the Legendre function, and the order r\ satisfies the equation 77(77+1) = 2 — ^ 
(Note that P v is the same for both roots of this equation.). We consider the case where 
>> 2, therefore rj « — | ± iy. Substituting the above expression for G m in Eq. (J3TJJ) 
gives 

I 1 h 2 G m ^(z \z(T))dT = A 2 I 1 fr 2 M {i*[cosh(As)]}dr . (D2) 

Here, s = s(z \z(r)). Changing the integration variable to s, and substituting <9 M P^ = 
s^d s (P^) in Eq. (jD2|) and integrating by parts gives 

J h 2 G m Jz \z(r))dr = -A 2 J d s (h 2 s,„-) P^cosh(Xs)]ds . (D3) 

Here s"i = s(zq\z{ji + e)). Note that both boundary terms vanished since h 2 {si) = 0, and the 
asymptotic form of P v (as s — > 00) [30| implies that asymptotically P'[cosh(As)] ~ e ^ 3 ^ 2 ^ Xs . 



We substitute 

.id 
ds 



P^[cosh(As)] = [Asinh(As)] -1 — P v 



21 



in Eq. (jD3|) . and integrate by parts once more. In the resultant equation we substitute the 
following expression for the Legendre function P v 



1 



P r? [cosh(As)] = — f= 



r(?7 + l/2) e 



(r?+l)As /111 i 



+ c.c. 



where F is the hypergeometric function. Noting that asymptotically 



r(z) 



r( z +i/2) 



~ z 



-1/2 



wc 



obtain 



lim 

m— >oc 



ri+e 



h 2 G m (zo\z(t))(1t = 



which conforms with our previous result. 
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